



Phenomenology of infrared finite gluon propagator and coupling constant 
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We report on some recent solutions of the Dyson-Schwinger equations for the infrared behavior of the gluon 
propagator and coupling constant, discussing their differences and proposing that these different behaviors can 
be tested through hadronic phenomenology. We discuss which kind of phenomenological tests can be applied to 
the gluon propagator and coupling constant, how sensitive they are to the infrared region of momenta and what 
specific solution is preferred by the experimental data. 
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I. INTRODUCTION 

During the last years there has been much effort in trying to 
obtain the infrared (IR) behavior of the Quantum Chromody- 
namics (QCD) Green's functions by means of theoretical and 
phenomenological studies CLEIMBB IaBl^l^llfllllll and 
by simulations of QCD on a lattice Qlofl^Tlal The in- 
frared behavior of the gluon and ghost propagators, as well as 
of the coupling constant are completely intertwined with the 
confinement problem. For example, if the gluon propagator 
would be as singular as 1 jk 4 when k 2 — ► 0, it would indicate 
an interquark potential rising linearly with the separation. 

In the late seventies, working in the Landau gauge and Eu- 
clidean space, Mandelstam obtained a solution of the Dyson- 
Schwinger equation (DSE) for the gluon propagator, in the 
case of pure gauge QCD, that behaved as 1/k 4 flol . A few 
years later Cornwall obtained a gauge invariant solution that 
behaved as 1 / '[k 2 + m 2 (k 2 )} 0. In this last case, as k — ► 0, 
the function m 2 (k 2 ) was interpreted as a dynamical gluon 
mass with the limit m 2 (k 2 — > 0) = m 2 . Both solutions repro- 
duce the expected perturbative behavior naturally at large k 2 . 
The 1 /k 4 result was named as a "confining solution", whereas 
the massive (or any infrared finite solution) became known as 
a "confined solution". It is clear that for the confined solution 
a more sophisticated explanation of quark confinement must 
be at work, and for this reason the appealing 1/k 4 solution 
became more popular, but not without some debate fl8il . 

Recently the coupled DSE equations for the gluon and 
ghost propagators were solved with different approximations 
121, resulting in an infrared gluon propagator that behaves 



roug hly as k 2 /[k 4 + a], i.e. vanishes as k 



0. 



This 

kind of behavior had also been predicted by Gribov 1 1 9tl and 
Zwanziger 1 20] in a different approach. The behavior of gluon 
and ghost propagators in Euclidean Yang-Mills theory quan- 
tized in the maximal Abelian gauge (MAG) were also stud- 
ied considering the effects arising from a treatment of Gri- 
bov copies in the MAG and those arising from a dynamical 
mass originating in a dimension two gluon condensate |21]. 
In these studies the infrared gluon propagator depends on the 
so called Gribov parameter, whose variation seems to make an 
interpolation between the vanishing propagator and the mas- 
sive one ! 

A great step ahead in this problem has also been pro- 
vided by the QCD lattice simulations in Landau gauge, which 
strongly support the existence of an infrared finite gluon prop- 



agator O2EIH0, where by finite we mean that the gluon 
propagator may be zero or different from zero at k 2 = 0. This 
is interesting enough, because indicates the appearance of a 
dynamical mass scale for the gluon, which imply in the ex- 
istence of a non-trivial QCD infrared fixed point, i.e. the 
freezing of the coupling constant at the origin of momenta 1 8 ] . 
Unfortunately the lattice data is not precise enough to defini- 
tively settle the questions of Green's functions at the origin of 
momenta, but there are claims that the results can nicely ac- 
commodate a massive gluon solution like the one proposed by 
Cornwall 11311 . Notice that at first glance it seems that these 
two possibilities are just a small detail, however they imply in 
two different confinement scenarios. Cornwall indicated that a 
dynamically generated gluon mass induces vortex solutions in 
the theory and these are responsible for the quark confinement 
IH7II . On the other hand the vanishing gluon propagator jointly 
with an infrared enhanced ghost propagator may induce an ef- 
fective linear confining potential between quarks 1 22] . It is not 
surprising that DSE for the QCD propagators lead to different 
solutions as long as they are solved with different truncations 
and approximations, as, for instance, the choice of the trilinear 
gluon vertex plays a crucial role in the solution. This means 
that we should wait for still more improved DSE solutions or 
lattice simulations for the gluon propagator. 



It is our purpose here to advocate that the hadronic phe- 
nomenology can provide solid information about the infrared 
behavior for the gluon propagator and coupling constant. 
Which kind of phenomenological tests are possible to propose 
in order to study the infrared behavior of the gluon propagator 
and coupling constant? In our opinion there are two possibil- 
ities: (i) We can make use of QCD inspired models where it 
is assumed "a priori" that the gluon propagator and running 
coupling constant must have one specific non-perturbative be- 
havior, as, for example, happens in the Pomeron model of 
Landshoff and Nachtmann 12311 : or (ii) we can consider per- 
turbative QCD calculations introducing an effective propaga- 
tor and coupling constant in the sense of the "Dynamical Per- 
turbation Theory" (DPT) proposed by Pagels and Stokar 
many years ago. In the next section we will detail and discuss 
examples of these two possibilities. 
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II. PHENOMENOLOGICAL TESTS FOR THE INFRARED 
BEHAVIOR OF THE GLUON PROPAGATOR AND 
COUPLING CONSTANT 

A. Testing the IR behavior in a QCD Pomeron model 

It is known that a connection between the theoretical pre- 
dictions of the perturbative Pomeron described by the BFKL 
Pomeron l25ll and the experimental data is far from being 
fully understood. Landshoff and Nachtmann proposed a sim- 
ple Pomeron model where the gluon propagation is affected 
by the non-trivial QCD vacuum leading naturally to a cor- 
relation length (or gluon mass scale). The existence of this 
correlation length is the only feasible explanation to the fact 
that diffractive interactions mediated by the Pomeron (which 
is assumed to be composed by at least two gluons in a color 
singlet state exchanged in channel t between hadrons) obey 
the additive quark rule. This model is quite suitable to test 
IR properties of coupling constant and gluon propagator ob- 
tained through DSE, because it always involves one integra- 
tion over the product of gluons propagators and coupling con- 
stant (g 2 (k 2 ) xD(k 2 )). 

As long as we deal with a QCD inspired model for the 
Pomeron we should be aware that the model will be valid up 
to a certain scale or approximation. Comparing the theoretical 
predictions of the Landshoff and Nachtmann (LN) Pomeron 
model with the experimental data we can be quite confident 
that the model is very reasonable to compute diffractive scat- 
tering in the small transferred momentum limit (i.e. small t 
physics). The Pomeron, in the LN model, is described as the 
singlet channel exchange of two non-perturbative gluons be- 
tween hadrons, and in the sequence we discuss the case of 
elastic differential cross-section for proton-proton scattering. 

In the LN model the elastic differential cross can be ob- 
tained from 



da 

dt 



16ns 2 



(1) 



where the amplitude for elastic proton-proton scattering via 
two-gluon exchange can be written as 



A(s,t) = is8a 2 [Ti-T 2 ] 



(2) 



with 



7i= jf d 2 kD(^+k)D(^-k)\G p {qM 2 (3) 



T 2 = [ d 2 kD(\+k)D(\-k) Gp ( q ,k-\ 



2G p (q,0)-Gp(q,k-- 



(4) 



where G p (q,k) is a convolution of proton wave functions 



Gp(q,k) = J d 2 pdK\\i*(K,p)\\i(K,p-k-Kq). (5) 



In Eq.0 D(q 2 ) is the non-perturbative expression for the 
gluon propagator. Note that the loop integral in Eq.0 
contains the product of two gluon propagators weighted by 
G p (q,k), which will be written in terms of proton form fac- 
tors. This is a general behavior, where the form factors will 
vary according to the type of hadron involved in the scattering. 
Evidently this means that the integral will depend on different 
scales contained in each form factor, which in the particular 
case of the proton is the proton mass, whereas for pion-proton, 
y — p scattering, or others, we have different mass scales in 
the integration, in addition to the one that enters in the gluon 
propagator and coupling constant. Another important feature 
of this model is the kinematic structure, where one of the glu- 
ons carry most of the momentum. The other gluon just seems 
to enter in the process with a very small momentum in order 
to form a color singlet, but it is responsible for the most im- 
portant part of the integration area 8261 . 

To obtain the cross section for the elastic differential p — p 
scattering we use the gluon propagator in Landau gauge writ- 
ten as 



(6) 



r 



where the expression for D(q 2 ) obtained by Cornwall is given 
by 



D- l ( q 2 )=[ q 2 +M 2 g ( q 2 )]bg 2 \n 
and the coupling constant equal to 



q 2 + 4M\ 



A 2 



ctsc(q 2 ) 



4k 



b\n[(q 2 +4Mj(q 2 ))/A 2 ] ' 



(7) 



(8) 



where M g (q 2 ) is a dynamical gluon mass given by, 



M\{q 2 )=m\ 



In 



q 2 +4m g 2 
A 2 



In 



-12/11 



(9) 



A(= Aqcd) is the QCD scale parameter. 
We will also make use of the running coupling constant ob- 
tained by Fischer and Alkofer 12711 which is given by 



a sA (x) 



oca(O) 



\n{e + a\x a i+bix b ^Y 



(10) 



where x — q 2 and 

0U(0) =2.972, «i =5.292 GeV~ 2fl2 , 



a 2 = 2.324, 



b\ =0.034 GeV 



-2bj 



b 2 = 3.169. and their respective propagator 
D(q 2 ) = Z(q 2 )/q 2 , where Z(q 2 ), in Landau gauge, is fitted 
by 



/™ I \ \ 1+28 



(11) 



2 



and 



R(x) 



dx 



2k 



1 



-dx 2K 



(12) 



k = 0.5953, 
c= 1.8934 GeV 



-2k 



where the constants appearing in Eg. I ll 1> and Eg. J12b are 

given by 

a^ 2 ) =0.9676, 
5 = -9/44, 
d = 4.6944 GeV~ 4K 

All these expressions were obtained from DSE solutions. 
Eqs.JlOb and were obtained in Landau gauge and corre- 
spond to a infrared vanishing gluon propagator. The gluon 
propagator comes with a definite mass scale once the renor- 
malization procedure is defined and the coupling constant 
fixed at one given scale. The Cornwall result was shown to be 
gauge invariant, and many of the features of this propagator 
and coupling constant are discussed by Aguilar and Papavas- 
siliou |28], its mass scale depends on the ratio m g /A. 

The details of the elastic differential p — p scattering calcu- 
lation can be found in Ref.|9], and the comparison of the re- 
sult with the experimental data of Breakstone et al. at s/s = 53 
GeV 1 29]is shown in Fig.Q. Note that the agreement between 
theory and experimental data is good only at small t, since at 
large t values we have contributions from 3-gluon exchange. 
We see in Fig.Q that Cornwall's propagator with a dynamical 
gluon mass of O(370)MeV fits the data quite well, whereas 
the Fischer and Alkofer one gives a quite large result for this 
cross section. We obtained the same result when analysing 
total hadron-hadron cross section |9, 11] and exclusive p pro- 
duction in deep inelastic scattering [9]. 

The integration of Eq.(|3} and Eq.© picks up an intermedi- 
ate region of momenta for the gluon propagators, i.e. above 
the region where the propagator that vanishes at the origin is 
suppressed and the coupling constant related to it is still large. 
As we shall see in the next examples most of the problems 
with this solution comes from the fact that the coupling con- 
stant given by Eq.(I10> is too large at the origin and up to sev- 
eral hundred MeV, and the integrations that we shall perform 
are not peaked at a momentum k 2 — 0, but are peaked in one 
region where the product g 2 (k 2 )D(k 2 ) appearing in Eq.Q and 
(0} is not small. As discussed in Ref.| 10] there are many in- 
dications that the experimental data is better described by a 
small coupling constant in the infrared. 



B. Testing the IR behavior in a perturbative QCD calculation - 
The y — > 7t° transition 

Perturbative QCD calculations involve the exchange of glu- 
ons at large momenta. These calculations do not provide very 
strong tests for the infrared behavior of the gluon propaga- 
tor, although a finite infrared gluon propagator gives a natu- 
ral cutoff for the many infrared divergent parton subprocess 
cross sections, which, in general, are dependent on the cutoff 
choice. On the other hand they may be influenced by the be- 
havior of the infrared coupling constant that comes out from 
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FIG. 1: Differential pp elastic cross section at v/(s) = 53 GeV com- 
puted within the Landshoff-Nachtmann model for the Pomeron, us- 
ing different infrared couplings and gluon propagators obtained from 
DSE solutions. 



DSE solutions, which at intermediate momenta may differ ap- 
preciably from the perturbative behavior. It is worth mention- 
ing that several calculations in the literature make use of in- 
frared finite coupling constants. For example, they are natu- 
ral in the so called Analytic Perturbation Theory (APT) l30ll 
or they appear as an effective charge in several perturbative 
QCD calculations l3~fll . In most of the cases these IR finite 
coupling are fundamental to confront the theoretical and ex- 
perimental data for several physical quantities. As mentioned 
before, we shall use the non-perturbative behavior of the cou- 
pling constant that comes out from DSE solutions in the sense 
prescribed by DPT Q. 

The photon-to-pion transition form factor F^Q 2 ) is mea- 
sured in single-tagged two-photon e + e~ — ► e + e~% Q reactions. 
The amplitude for this process has the factorized form 



FyxiQ 2 ) = -1^ dx^ K (x,Q 2 )T^(x,Q 2 ), (13) 
where the hard scattering amplitude T^.(x,Q 2 ) is given by 



(14) 



Using an asymptotic form for the pion distribution amplitude 
<Ptu = V3fnx(l —x), we obtain 1 10] 



Q%(G 2 ) = 2A( i-I^P 



(15) 



where Q* = exp 3 / 2 Q is the estimated Brodsky-Lepage-Mac- 
kenzie scale for the pion form factor in the scheme discussed 
inRef.il. 
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In Fig.© we compare the photon to pion transition form 
factor with CLEO data 1 33 ] . The curves were computed with 
different expressions for the infrared behavior of the running 
coupling constant. We assumed f % ~ 93 MeV and A = 300 
MeV. We also made use of a running coupling constant deter- 
mined by Bloch 1 34], where the absence of the Landau pole at 
q 2 = A 2 is reminiscent of APT, which is given by 



a sB {q 2 ) 



1 



(16) 



co - 



r- 



CQOCo + ■ 



471 



Po Vi°g(0 *-i 



where / = g 2 /Ag C£) , co = 15, (Xo = 2.6, and po = H — 
where rif is the number of flavors. Ea.( ll6> is also consistent 
with a propagator that vanishes at k 2 = 0. With the coupling 
constant of Ea. (ll6> we obtain a fit for the photon-pion transi- 
tion form factor very far from the experimental data. The re- 
sult obtained when we use Eq.dlOl is not shown and gives an 
even worse fit. The infrared value of the coupling constant is 
so large in the case of the coupling constants given by Eqs.( I10> 
and ( I16> . that we are not sure that the perturbative result can 
be trusted even at such large momentum scale. The momen- 
tum scale appearing in Fig.(0 is above the GeV scale, where 
the DSE solutions for the gluon propagators already assume 
their perturbative 1/k 2 behavior. The values of the infrared 
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expressions for the infrared behavior of the running coupling con- 
stant. 



coupling constants related to the class of DSE solutions con- 
sistent with a vanishing gluon propagator are much stronger 
than most of the phenomenological estimates of the frozen 
a,(0) value that we quoted in Ref.|10] (a,(0) w 0.7 ±0.3), 
and are at the origin of the strange lower curve of Fig.0. The 
data is only compatible with Eq.(|8}, which has a smoother in- 
crease towards the infrared region. Perhaps this behavior is 
actually indicating that the transition to the infrared should be 
a soft one. 

Note that in Fig. 10 the curves obtained with Cornwall's 
coupling constant do not show large variation in the full range 



of uncertainty of the dynamical gluon mass. It is interesting 
that its behavior is quite stable in this case as well as for the 
pion form factor studied in Ref. llOll . If we had large varia- 
tions of the infrared coupling constant with the gluon mass 
scale we could hardly propose any reliable phenomenological 
test for its freezing value. We stress that the results are ob- 
tained for a perturbative scale of momenta and it seems that 
we have to choose between two possibilities: (f) This process 
cannot be predicted by perturbative QCD up to a scale of sev- 
eral GeV, or (if) some of the DSE solutions are predicting a 
too large value of the coupling constant in the infrared and 
the approximations made to determine these solutions are too 
crude. 

We have also computed the perturbative pion form factor 
(F x (q 2 )), and the same happens in that case, i.e. only the so- 
lutions for the gluon propagator and coupling constant deter- 
mined by Cornwall, when plugged into the perturbative QCD 
expression for F K , produce the match between the theory and 
the experimental data. The details of the pion-photon tran- 
sition form factor calculation can be seen in Ref.|9] and the 
pion form factor calculation is described in Ref.[ 10]. 



C. Testing the IR behavior in a perturbative QCD calculation - 
Total hadronic cross sections within the parton model 

Total hadronic and jet cross sections can be calculated in a 
straightforward way within the parton model (see Ref.[35] for 
a review); The cross section for producing jets with pj > pT„, in 
through the dominant process gg — ► gg is given by 



dp 2 ddgs 



dpj 

dx 1 dx 2 g(xuQ 2 )g(x2,Q 2 ), (17) 



X\Xi>Apf/s 

where g(x,Q 2 ) is the gluon flux and p\ is the momentum 
above which we can use the perturbative calculation of the 
subprocess differential cross section da gg /dpj. This momen- 
tum scale is usually assumed to be larger than 1 GeV 2 . At 
high energies, the first order perturbative result for the cross 
section a gg can be written in terms of the variable § as ]|3F 



9kOq 



(18) 



where mo is a particle production threshold and oco is an effec- 
tive value of the running coupling constant. This cross section 
causes a rapid increase in Oj et (s) if a eikonalization procedure 
is not used in the calculation. Assuming the simple ansatz for 
the gluon flux 



(I-*) 5 



a straightforward calculation yields 



'jet 



97CCX0 



(19) 



(20) 
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where e = J — 1 > 0. In this calculation we have neglected 
some factors in the right hand side of the final expression. 
Nevertheless, in the limit of large enough s this expression 
reproduces the expected asymptotic energy dependence of 
Gjet(s). Moreover, with specific values for mo and (Xo, it is 
possible to show that at ^/s ~ 630 GeV, this jet cross section 
is of order of the pp and pp total cross sections l36ll . It is also 
clear that Eq. J20> does violate unitarity. 

The simple derivation of the cross section behavior that we 
have seen above is not too useful due to the following rea- 
sons: a) the m\ and the terms in Eq. \20l are totally ad hoc 
and b) Eq. J20> violates unitarity. Unitarity is recovered with 
the eikonalization of the model. However the procedure still 
keeps the ad hoc constants Wq and 0Cq as parameters that can 
be obtained only with the data fitting. In Ref. II37II the ele- 
mentary gluon-gluon cross section was calculated within the 
dynamical perturbation theory scheme (DPT) 1 24], where the 
effective gluon propagator and coupling constant enters into 
the calculation. With this procedure we eliminate the freedom 
existent in the previous calculation, due to the choice of Mq 
and 0Cq, changing hjq by m 2 and still have a parameter less 
in the calculation because the running coupling now is also a 
function of m 2 . 

In the eikonal representation the total cross section is given 



a tot (s) =4ti / dbb[l-e-*'^cos% R (b,s)}, (21) 



by 



where s is the square of the total center-of-mass energy, b 
is the impact parameter, and %(b,s) = % R (b,s) +i%j(b,s) is 
a complex eikonal function. In the QCD eikonal model with a 
dynamical gluon mass, henceforth referred to as DGM model, 
we write the even eikonal as the sum of gluon-gluon, quark- 
gluon, and quark-quark contributions: 



%qq(b,s) and %q g (b,s), needed to describe the low-energy for- 
ward data, are parametrized with terms dictated by the Regge 
phenomenology: 



X qq (b,s) = i-LA-^W(b;n 9g ), 



(25) 




W(b;^JI^), (26) 



where A, A', B' , n qq and p. gg are fitting parameters. 

The innovation in our approach is that the subprocesses 
cross sections that appear in Ea.J22> are computed within 
DPT, where the gluon propagator and coupling constants are 
the ones obtained from DSE solutions. For example, O gg (s) is 
computed with the help of the Cornwall propagator and cou- 
pling constant, where, for simplicity, we neglected the mo- 
mentum dependence of the running mass in the denominator 
of the gluon propagator, obtaining 



3%a 2 . 



-31n 



12.f 



s4 



- 55m 2 s 3 + 12m 4 g f + 66mp- 



4tn 2 s[s — m 2 ] 2 



3mi 



m 2 



(27) 



.25 r 



O 
Q 



•2 - • 



% + (b,s) = % qq (b,s)+% qg (b,s)+% gg (b,s) 

= i[a qq {s)W{b;n qq )+a qg {s)W{b;n qg ) 
+ o gg (s)W(b;fi gg )], 



(22) 



1.15 



where % l P p(b,s) = % + (b,s) ±%~ (b,s). Here W{b;^i) is the 
overlap function in the impact parameter space and Cij(s) are 
the elementary subprocess cross sections of colliding quarks 
and gluons = q,g). The overlap function is associated 
with the Fourier transform of a dipole form factor, 



W ( b 'V) = -^ l {vbfK,{pb) 1 



(23) 



where Kj (x) is the modified Bessel function of second kind. 
The odd eikonal %~(b,s), that accounts for the difference be- 
tween pp and pp channels, is parametrized as 



X-(b,s)=C-Z T ^e^ 4 W(b;pr), 



(24) 



where m g is the dynamical gluon mass and the parameters 
C and /j~ are constants to be fitted. The eikonal functions 



1.1 
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FIG. 3: The % 2 /DOF as a function of dynamical gluon mass m g . 

The full detail of the calculation can be found in Ref. 
We fitted the pp and pp scattering data keeping m g as a 
free parameter. Our global fit results indicate a minimum 
value just about m g m 400 MeV. These results are shown in 
Fig-@> where a general dashed curve is added to guide the 
eye. Roughly, taking a 5% variation on the minimal % 2 /DOF 
value indicated by the general curve, it is possible to estimate 
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a dynamical gluon mass m g w 400_jqq MeV. This value ob- 
tained through data fitting is totally compatible with the ones 
found by Cornwall and subsequent determinations ma 
The fits for the total cross sections (o to t) can be seen in Fig.(0} 
in the case of m g = 400 MeV. 
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FIG. 4: Total cross section for pp (solid curve) and pp (dashed curve) 
scattering. 

It must be noticed that in the cross sections calculations we 
do have another integration over the gluon propagators and 
coupling constant expressions weighted by the parton struc- 
ture functions. In this case we also test the infrared behavior 
of the gluon propagator because a large part of high energy be- 
havior of the total hadronic cross section is due to soft gluons 
l35ll . which means small £ values in Ea. (l27> corresponding 
to the IR contribution of the gluon propagator. An analysis 
similar to the pp total cross section calculation shown is this 
subsection has also been performed in the case of y — p and 
y — y scattering with the same dynamical gluon mass 1 38]. 



are originated from different approximations or truncations of 
the DSE. As one example we could recall that for many years 
an IR DSE solution for the gluon propagator behaving like 
1 /k 4 was very popular, but recent lattice QCD simulations ba- 
sically discarded such possibility. 

Lattice QCD simulations have already proportioned a great 
improvement in our knowledgement of infrared Green's func- 
tions. There are strong indications that the gluon propagator 
is finite and possibly the same may happens for the coupling 
constant. This result is important because it may give some in- 
formation about the confinement mechanism. Unfortunately 
the lattice result is still not precise enough near the origin of 
momenta in order to discriminate between the different gluon 
propagator solutions. 

In this work we are advocating that hadronic phenomenol- 
ogy can distinguish between the different IR behaviors of 
DSE. No matter how we deal with QCD inspired models, like 
the LN Pomeron model, or just perturbative QCD calculations 
improved by the use of an effective gluon propagator or cou- 
pling constant, we see that Cornwall's propagator is selected 
by the experimental data. In all the examples discussed in the 
previous section, we see that the phenomenological informa- 
tion is non-trivial, in the sense that it results from the calcula- 
tion of physical quantities where the gluon propagator or prod- 
uct of propagators are integrated weighted by different func- 
tions (involving different mass scales), and all quantities show 
agreement with the experimental data for gluon masses that 
are in the same range of masses predicted by Cornwall several 
years ago. It is hard to believe that such coincidence is a for- 
tuitous one. These values for the gluon mass lead to a frozen 
coupling constant at the origin of momenta whose value is 
not larger than O(l). Certainly an infrared finite coupling 
constant is welcome for hadronic phenomenology 
Maybe future lattice calculations will shed some light on this 
problem, indicating which DSE approximation is reasonable 
or not, which confinement scenario is more probable to be at 
work, and which kind of improved phenomenological calcu- 
lation can be done with such infrared finite gluon propagator 
and coupling constant. 



III. CONCLUSIONS 

We discussed DSE solutions for the gluon propagator. 
These solutions can be divided in two possibilities according 
to their behavior as k 2 — ► 0, one vanishes at the origin and 
the other has a value different from zero and consistent with 
a dynamically generated gluon mass. We should not be sur- 
prised with the existence of different solutions because they 
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